Fractional derivatives in Dengue epidemics* 
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Abstract 

We introduce the use of fractional calculus, i.e., the use of integrals and derivatives of 
m : non-integer (arbitrary) order, in epidemiology. The proposed approach is illustrated with an 

outbreak of dengue disease, which is motivated by the first dengue epidemic ever recorded in 
the Cape Verde islands off the coast of west Africa, in 2009. Numerical simulations show that 
in some cases the fractional models fit better the reality when compared with the standard 
differential models. The classical results are obtained as particular cases by considering the 
order of the derivatives to take an integer value. 

Keywords: fractional calculus, Riemann-Liouville derivatives, epidemiological model, dengue 
disease, numerical approximation. 
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1 Introduction 

The incidence of dengue has grown intensely in the last decades and, according to the World Health 
Organization, about 40% of world's population is now at risk. This global pandemic is attributed to 
the unprecedented population growth, the rising level of urbanization without adequate domestic 
water supplies, increasing movement of the virus between humans (due to tourism, migration, or 
international trade), and lack of effective mosquito control pQ. Dengue virus is transmitted to 
humans through the bite of infected Aedes mosquitoes, specially Aedes Aegypti. Once infected, a 
mosquito remains infected for life, transmitting the virus to susceptible individuals during feed. 
Without a vaccine, vector control remains the only available strategy against dengue. 

Appropriate mathematical models can give a deeper insight into the mechanism of disease 
transmission 2, 3, 4 . Typically, epidemiologic models are formulated with classical derivatives of 
integer order. Here we propose the use of generalized fractional derivatives. Fractional calculus 
(calculus of non-integer order), in spite of its long history as a pure branch of mathematics, only 
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recently has shown to be useful as a practical tool [5]. In this paper we claim that fractional 
calculus provides an interesting modeling technique in the context of epidemiology. 

We begin by considering a simple epidemiological model that represents an episode of dengue 
disease. The rest of the paper is then dedicated to introduce the notion of fractional derivative 
in the sense of Riemann-Liouville, reformulating the dynamics of the classical model in terms 
of fractional derivatives, and finally applying a recent approximate technique to obtain numer- 
ical solutions to the fractional model. Numerical simulations show that models with fractional 
derivatives may provide a better description of reality when compared with the standard ones. 

2 Epidemiological model 

In this paper we assume that the host population is divided into three classes: susceptible, Sh(t), 
individuals who can contract the disease; infected, Ih(t), individuals capable of transmitting the 
disease to others; and resistant, Rh(t), individuals who have acquired immunity at time t. The 
total number of hosts is constant, i.e., Nh = Sh(t) + Ih{t) + Rh{t)- Similarly, we also have two 
compartments for the mosquito: S m (t) and I rn (t) with N m — S m (t) + I m (t). The model is 
described by the system of differential equations 



subject to given initial conditions S/j(0), Ifc(O), Rh(0), S m (0) and 7 m (0). The recruitment rate of 
human and vector populations are denoted as [ihNh and /i m iV m , respectively. The natural death 
rate for humans and mosquitoes is described by the parameters /z./, and /x m , respectively. We 
assume that B is the average daily biting (per day) of the mosquito whereas j3 m h and fthm are 
related to the transmission probability (per bite) from infected mosquitoes to humans and vice 
versa. The recovery rate of the human population is denoted by r\h ■ 

3 Fractional Calculus 

Fractional calculus was originated at the end of the seventeenth century and consists in the study 
of derivatives and integrals of arbitrarily real or even complex order. Recently, fractional calculus 
is experiencing an intensive progress in both theory and applications [3] . Main claim is that 
a fractional model can give a more realistic interpretation of natural phenomena. Among several 
different definitions that can be found in the literature for fractional derivative, one of the most 
popular is the Riemann-Liouville derivative. 

Definition 1 Let x(-) be an absolutely continuous function in [a,b] and < a < 1. The (left) 
Riemann-Liouville fractional derivative of order a, a D^ , is given by 



We reformulate system (flj using fractional derivatives. To this end, we simply substitute the 
first-order derivatives by Riemann-Liouville derivatives of order a: 
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Note that when a — y 1 the fractional system @ reduces to (JXJ) . Given a certain reality, our goal 
is to find the order a that makes the model more realistic. As we shall see in the next section, the 
best value of a in system @ is in general different from one, i.e., the classical model ([T]) is often 
not the best choice. 



4 Numerical simulations 

Similarly to the classical theory of differential equations, there are no general methods to solve 
systems of fractional differential equations analytically. Furthermore, the fractional case is be- 
lieved to be more difficult to handle even numerically [7^ • Recently, an approximation based on a 
continuous expansion formula for the Riemann-Liouville fractional derivative has been proposed 
in [8] and further improved in [9]. The basic idea is to express fractional terms by means of a 
series involving integer-order derivatives. The immediate result is that one can then use classical 
methods to obtain an approximate solution to the original fractional problem. With the same 
assumptions of Definition [I] the approximation for a Df x(-) is 

N 

a D?x(t) ~ A(a,N)(t - a)- a x{t) + A'(a,N)(t - a) 1 "^) -^C(a,p)(i - a) 1 -*- a V p (t), (3) 

p=2 

where 

V p (a) = 0, p = 2,3,...,N, 
and A = A(a,N), A' = A'(a,N) and C p = C(a,p) are given by 
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Using ([3]), we approximate system @ by a system of ordinary differential equations in which, 
subject to the order of approximation N, new variables of the form Vp(-) will appear: 



(t) = \vi h N h - (B(3 mh j^ + fi h )S h - At- a S h + Ep= 2 C p t x -P- a V p s ^(t) 

f^(t) = (1 - p)(t - a)P- 2 S h (t), p = 2, 3, . . . , N 
%-{t) = \Bp mh ^-S h - ( Vh + fj, h )I h - At- a I h + E^ =2 C^-v-^it) 



dS h 
dt 

dV T ' 
dt 



A'-H c 



A'~H C 



dV, 



dt 
dR 



£-(*) = (1 - P )(t - a)P- 2 I h (t), p = 2, 3, . . . , N 



^t(t) = r) h I h - fi h R h - At- a R h + Ep =2 Cp^-P-oV** (t) 
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A'~H 



dt 
dV, 



(t) = (1 - p)(t - a)P~ 2 R h (t), p = 2, 3, . . . , N 



V m N m - (Bp hm jft + ^ m )S m - At- a S m + Ep =2 Cpti-P-^Vp 3 ™ (t) 



dt 



(t) = (1 - p)(t - a)P- 2 S m (t) 7 p = 2, 3, . . . , N 



dl,„ 
dt 
dV' 



(*) = Bp hm ^-s m - (fi m )i m - At- a i m + j:^ =2 c p t l -p- a v p ^ (t) 



A'-H - 1 



fe-(t) = (1 - p)(t - a)P- 2 I m (t), p = 2, 3, . . . , N. 



dt 



(4) 

System ([4J is solved numerically using the following values that are inspired in the outbreak 
occurred in 2009 in Cape Verde: N h = 56000, B = 0.7, (3 mh = 0.36, /3 hm = 0.36, /j, h = 1/(71 x 365), 
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Vh = 1/3, fJ, m = 1/10, m — 3, and N m — my. iV/,. The initial conditions for the system ([l} of 
ordinary differential equations are: 5^(0) = Nh — 216, lh(0) — 216, Rh(0) = 0, 5 m (0) = N m and 
^m(O) = 0. These values, together with a zero initial condition for Vp h , Vj h , V p Rh , Vp m , and V p Im , 
make system @ to be an ordinary initial value problem that can be treated numerically. Figure [TJ 
shows the solution to the systems (JTJ and (QJ with respect to the variable The system (QJ is 
solved with N = 7 and different values of a. In order to compare both fractional (a G [0, 1)) and 




Figure 1: solution to the classical model (JTJ) ver- Figure 2: real data versus the solution to the 
sus the solution to the fractional model © with classical model (JTJ and the solution to fractional 
different values of a. model (JSJ with a — 0.987. 

classical (a = 1) models, we include the statistical information from Cape Verde in Figure O The 
best model is, in our case, the fractional model with a = 0.987, which corresponds a percentage 
error of thirteen. In contrast, the percentage error associated with the classical model is sixty two. 



5 Conclusions 

Describing the reality through a mathematical model, usually a system of differential equations, 
is a hard task that has an inherent compromise between simplicity and accuracy. In this paper we 
consider a very basic model to dengue epidemics. It turns out that, in general, this basic/classical 
model does not provide enough good results. In order to have better results, that fit the reality, 
more specific and complicated set of differential equations have been investigated in the literature 
- see [TOJ [13 [TJ] and references therein. Here we propose a completely new approach to the 
subject. We keep the simple model and substitute the usual (local) derivatives by (nonlocal) 
fractional differentiation. The use of fractional derivatives allow us to model memory effects, and 
result in a more powerful approach to epidemiological models: one can then design the order a of 
fractional differentiation that best corresponds to reality. The classical case is recovered by taking 
the limit when a goes to one. Our investigations show that even a simple fractional model may 
give surprisingly good results. However, the transformation of a classical model into a fractional 
one makes it very sensitive to the order of differentiation a: a small change in a may result in a 
big change in the final result. 

The present work can be extended in several ways: by fractionalizing more sophisticated 
models; by considering different orders of fractional derivatives for each one of the state variables, 
i.e., models of non-commensurate order. 
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